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MACROSCOPIC FLUCTUATIONS THEORY OF AEROGEL 

DYNAMICS 
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Abstract. We consider the thermodynamic potential describing the macroscopic 
fluctuation of the current and local energy of a general class of Hamiltonian models 
£^ ■ including aerogels. We argue that this potential is neither analytic nor strictly 

convex, a property that should be expected in general but missing from models 
studied in the literature. This opens the possibility of describing in terms of a 
thermodynamic potential non-equilibrium phase transitions in a concrete physical 
context. This special behaviour of the thermodynamic potential is caused by the 
fact that the energy current is carried by particles which may have arbitrary low 
'-"H . speed with sufficiently large probability. 

1. Introduction 

^-j- ! In equilibrium statistical mechanics, the Boltzmann-Gibbs probability distribu- 

tion is at the heart of the connection between the micro- and macroscopic level of 
systems made of a large number of particles N. From that distribution, one may 
compute the relevant thermodynamic potentials which may be interpreted as rate 
functions at speed N governing the large deviations of the macroscopic quantity of 
interest. Typically, the lack of strict convexity of the thermodynamic potentials is 
the signature of the existence of phase transitions. Out of equilibrium, there is not a 
such thing as the Boltzmann-Gibbs formula for the microscopic configurations and, 
in order to understand the macroscopic behaviour of a given material, one must pro- 
ceed differently. Building on the seminal work of Onsager and Machlup [23] , several 
d ' authors have built a general picture of the macroscopic current fluctuations in non- 
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equilibrium statistical mechanics, see [31 [6] for overviews of the subject. Although 
this approach is rigorously justified only in the framework of stochastic interacting 
particles, the scope is much more ambitious. Macroscopic fluctuation theory de- 
scribes, after a diffusive rescaling of space and time, the macroscopic behaviour of 
extended systems crossed by a flow of a conserved quantity: energy, mass or other 
relevant quantities. In this approach, the classical thermo dynamical potentials are 
replaced by a rate function describing the fluctuations in space-time of the local cur- 
rents and energies. From the point of view of the theory of large deviations this rate 
functional is strictly analogous to the equilibrium thermodynamic potentials in the 
sense that it is a large deviations rate function at speed N. From the knowledge of 
this rate function one may, as in equilibrium, compute the cumulants of the relevant 
physical quantities. Starting from a quadratic form for the rate functionnal, this is 
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what is achieved in [5] where a relation for the successive cumulants of the current is 
obtained. A fundamental issue, on which almost no result exist for realistic physical 
systems, is to determine the form of the non-equilibrium rate function. 

In this paper, our goal is to derive the form that this rate function must take in 
a general class of Hamiltonian models including aerogels. These materials are gels 
whose liquid component has been removed and replaced by atoms of gases. Because 
of their particular insulating properties, they have proved useful in applications 
ranging from the manufacturing of new insulating glasses to the engineering of key 
parts of space probes. The theoretical models we will consider are easy to study 
numerically and well described by a local equilibrium assumption [TU [15] in a weakly 
interacting regime. Using a stochastic approximation, we show that the function 
describing the macroscopic fluctuations of those models is drastically different from 
the one generally expected and obtained since the work of Onsager and Machlup [TJ 
[21 El IH El [71 [23]. As a function of the current it is not a quadratic function around 
the stationary current. More than that, it is neither a smooth nor a strictly convex 
function. This means, by analogy with the equilibrium case, that aerogels must 
exhibit interesting out of equilibrium phase transitions phenomena. This should 
also apply to a class of models made of tracers and scatterers recently considered 
by several authors [121 EH El 1201 [21] . 

We emphasize that our main working hypothesis is the validity of the descrip- 
tion of the systems at a diffusive space-time scale. This amounts to assuming the 
existence of a rate function and does not have any implication a priori on the spe- 
cific shape that this function must take. The validity of the diffusive description of 
macroscopic systems should hold for a large class of mechanical systems. On the 
other hand, there are very few a priori constraints imposed on the rate function 
except a global one imposed by the Gallavotti-Cohen symmetry and arguably local 
ones (on the first two derivatives) imposed by linear response theory. We will see, 
however, that the collisional mechanism by which the energy is transported through 
the system by the microscopic components is responsible for a singular shape of the 
rate function. The two constraints stated above are still satisfied. 

2. The models 

The dynamics that we want to consider is as follows. Consider N particles of unit 
mass with positions and momenta (q, p) = {(qj, Pi)} 1< . <JV , with qj,Pi E M d . The 
positions are measured with respect to N fixed centers located on a ID lattice. The 
Hamiltonian H takes the form 



N 



ff(p,s) = £ 



=i l 



^ + V(q,) + U(q t - q i+1 ] 



(2.1) 



where the interaction potential U is equal to zero inside a region Qu C Mr with 
smooth boundary A of dimension d — 1, and equal to infinity outside. Likewise, 
the pinning potential V is assumed to be zero inside a bounded region fly an d 
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infinite outside, implying that the motion of a single particle remains confined for 
all times. The regions Qjj and Qy being specified, the dynamics is equivalent to a 
billiard in high dimension. A typical example is given by the figure below. Each 
disk moves freely within its square cell and collides with a neighbor when they both 
get sufficiently close to the hole located in the wall separating their adjacent cells. 









Figure 1. Simplified aerogel dynamics 

The first type of such models, in which the local dynamics is described by semi- 
dispersing billiards, was introduced in [10] and its thermal transport properties were 
studied in [13]. General collisional models of this type have been introduced and 
analyzed in [15]. To model those systems out of equilibrium, one couples them at 
their boundaries with stochastic thermal walls at different temperatures. Rigorously, 
almost nothing is known about the invariant state of such systems. Heuristically, 
however, the behavior of such systems is easy to guess [TSJ [TS]. The evolution of the 
energy of the n-th particle may be written as, 

E n (t) - E n (0) = J n _i^ n [0, t] - J n ^ n+ i[0, t}. (2.2) 

where 

J n _ n+1 [0,t] = \Y, [Pn(s k ) 2 -p^ +1 (sk) 2 }, (2.3) 

0<k<N t 

and p^ is the component of the vector p n in the direction of the unit vector fi = 
\\qi — gj + i|| _1 (gj — qi+i) at the time of collision, namely p^ = p n • n. The integer 
N t counts the number of collisions up to time t and Sk is the k-th collision time. 
Assume now that such systems are thermalized at different temperatures at their 
boundaries. In order to understand the transfer of energy from one side to the 
other, one is interested in the ergodic behavior of the current and in computing 
lim t ^ +00 1 _1 J n ^ n+ i[0, t], giving the average current of energy in the stationary state. 
Because of the special form of the time- integrated current ( 12. 3p . a natural guess to 
make is that local equilibrium settles in and to assume that this limit is given by 

lim -J n ^ n+1 [0,t] = v(T n -T n+l ) (2.4) 

where v = lim^oo £ _1 iV t is the frequency of collisions between neighbors under local 
equilibrium conditions and T n = |(pj[) is the average kinetic energy of the particles. 
The conductivity is thus identified to the frequency of collisions. Typically, the 
collisions occur when the particles get near the boundaries of their cell and thus 
the frequency of collisions is roughly proportional to the average number of visits 



4 RAPHAEL LEFEVERE, MAURO MARIANI, AND LORENZO ZAMBOTTI 

to the boundaries per unit time. Because the particle travels freely within its cell, 
this is proportional to the square root of the average energy and thus to y/T^. 
Numerical studies show that the identification of the conductivity with the frequency 
of collisions holds true to a very high degree of accuracy in a wide class of collisional 
dynamics, when the individual particles collide rarely [T4"l IT5]. 

3. Macroscopic fluctuation theory 

In macroscopic fluctuation theory, the evolution of the relevant variables is de- 
scribed after a diffusive rescaling. One introduces a macroscopic evolution time 
variable s = t/N 2 and a continuous space variable x to describe functions on the 
spatial lattice by functions on the interval [0, 1]. One assumes that the system may 
be described by smooth functions J and e related to the original time-integrated 
current and energy by 

■Wi[o,t] = ffj{jj,^ 

We define the instantaneous current at the macroscopic scale by 

j(x,s) = d s J(x,s) 

and we introduce the notation: jk(s) : = j (77, <s) • The evolution of the energy is 
governed by (12.21) which becomes at the macroscopic scale: 

d s e(x, s) = -d x j(x, s). (3.2) 

Macroscopic fluctuation theory predicts that asymptotically in N for j and e related 
bydS2D 

F({e(x,s'),j(x,s')}) ~exp[-iVX(j,e)] (3.3) 

where X(j, e) is given by 

X(j,e) = [ ds' [ dxct>{j{x,s'),d x e{x,s'),e{x,s')) (3.4) 

Jo Jo 

if j and e satisfy (13.21) . and X = +00 otherwise. The shape of is a priori quite 
arbitrary except for the constraints mentioned in the introduction. Nevertheless, in 
the stochastic models that have been studied so far, it is known that is a quadratic 
function of the current around the stationary current in the diffusive scaling limit 
[H El El El El E]- We assume that macroscopic fluctuation theory applies to our 
models in the sense that there exist smooth functions e and J and that (13. 3p and 
(13. 4p hold true. We argue now that in the deterministic models described above <fi 
can not be quadratic in j and is rather given by (13.71) below. We assume that the 
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functions <fr, j and e are sufficiently smooth, we discretize space with increments of 
size 1/N and time with small steps of size As and set: 

a [a/As] N 

Z(j,e) = lim lim — J^ J2<t>k,i (3.5) 

Z=l fe=l 

With, 

M := (j) (j k (lAs), N[E k+l (lAs) - E k (lAs)\, E k (lAs)) 

and we identify exp{— As^;} as the probability of observing an instantaneous 
current j k (lAs) at site k during a small macroscopic time interval As, while the 
local energy and energy gradient are fixed. A key observation for our purpose is 
the following: as the macroscopic time varies by a small increment As, the energy 
e changes only by an amount proportional to that increment. On the other hand, 
in terms of the original microscopic variables, this corresponds to very long time 
scale iV 2 As. Observing an instantaneous current j k {s) over a small interval of time 
As corresponds to observing an average time-integrated current J k [N 2 s, N 2 (s + 
As)] = jk{s)NAs. Because of this relation we identify the probability of finding an 
instantaneous current j k (lAs) with the probability of observing an average time- 
integrated current (over a microscopic time interval iV 2 As) equal to jk{lAs)/N. 
Thus, 

*W 1 , „/.4[iV 2 (As,iV 2 ((+l)A S )] j k (lAs) 



tp ^A S l0gP ( — W. --"it) (3 - 6) 

For iV large, the right-hand side converges to the large deviation functional of the 
time-integrated current J evaluated at jk/N, with the energy variables fixed. In our 
context, we are led to the following picture: a given particle makes a huge number 
of collisions with its neighbors while the average energy of its neighbors is basically 
unaffected. We assume further that the dynamics is sufficiently chaotic so that 
the dynamical correlations between neighbors are destroyed between two collisions 
among them. This leads us to model the system during a time slice As by the 
following stochastic dynamics. Each particle moves freely in its cell and the effect 
of the collisions of a given particle with its neighbors is modeled by the action of 
hot walls at fixed temperatures. The temperatures of the walls correspond to the 
average energies of the particles, which are fixed during the small macroscopic time 
interval As. Pictorially, the model of figure 1 is transformed into the model of the 
figure below. 



^ ■ ^H ^ ■ ^ /J 



Figure 2. Collisions replaced by interactions with hot walls 
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Clearly, the direction orthogonal to the array of heat baths is irrelevant to the 
transport of energy and one is thus lead to the model of confined tracers introduced 
in [19]. Each particle moves freely in an interval and when it encounters the hot 
walls located at its boundaries, its velocity is reflected and it gets a new random 
speed v distributed according to the distribution: 

2 

(pp(v) = (3v e - ^, v > 0, 

where (3 is the inverse of the temperature of the wall where the collision has taken 
place. We note that when the temperatures on the boundaries are equal, the in- 
variant distribution of the random process describing the evolution of the position 
and velocity of the particle is given by the equilibrium Maxwell-Boltzmann-Gibbs 
formula [12] . It is very important to note that if one takes as the updating rule a 
Maxwellian distribution, one does not get the Maxwell-Boltzmann-Gibbs distribu- 
tion as the invariant measure for the random process □• In the stochastic model the 
energy current corresponding to a given particle is given by: 

1 Nt 
J[0,t] = a -J2(-l) k v 2 k 

k=\ 

where Vk is the speed the particle gets at the fc-th collision with a wall, N t is the 
number of collisions up to time t and <7o is the sign of the initial velocity. It is impor- 
tant to note at this stage that the form of the current is based on the deterministic 
part of dynamics, it simply corresponds to the kinetic energy carried by particles 
moving back and forth. As we will see below, Nt, whose large deviation properties 
will be central to our argument, depends in a simple fashion on the sequence of 
velocities assumed by the particle. This is because the deterministic motion of the 
particle between the interactions with the walls is ballistic. One can show that, 

Um /[<M]_ T L -T R 



\2T L ' ~ V 2T R ' 



where Tjr, and Tr are the left and right temperatures. 

In [18] , we study in details the large deviations functional of the current I(j, r, T) 
defined as: 

■^[0'*] „\ „-tl(j,r,T) 



P Tj r is the probability associated to the stochastic dynamics with a fixed temperature 
difference r = Tl — Tr between the two walls and average temperature T = L ~t R . 



In the context of the study of gases contained in a vessel, it is customary to model the action 
of thermal walls by imposing special boundary conditions. We show in [18] that this amounts to 
specify the domain of the generator of the random process corresponding to the above updating 
rule. 
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More precisely, we show in [T8] the following result: if r 7^ then, 



\ime- 2 l(ej,ET,T) 



-j^- if JT>KT 2 
if JT 6 [0, KT 2 } 

-& if jte[-kt 2 ,0] (3.7) 



I ^f- if JT < -KT 2 

where k = (|j;) 5 - When r = 0, hm e ^ £ ~ 2 %( £ j, e ^i T) = -±tfz- Denoting the RHS by 
G(J,t,T), we observe that it has the Gallavotti-Cohen symmetry: 

g^ r ,T)-g(-j,T,T) 



2T 2- 

The most striking feature of the functionals X(-, r, T) and Q(-,r,T) is the fact that 
they vanish on the interval [0, kt]. The origin of this phenomenon may be traced 
back to the occurrence of slow velocities with sufficiently large probability in the 
course of the motion of the particle. More precisely, the process can be entirely 
characterized by the sequence (tfc)fceN of waiting times between collisions of the 
particle with the walls. Those waiting times are equal to the inverse of the random 
speeds t k = l/u* and are thus distributed according to a density: 

(3 _ n 

The study of the large deviations of the current boils down to analyzing the large 
deviations of the Markov renewal process N t = sup{n : S n < t}, where S = and 
S n = ti + ... + t n . The fact that the distribution i^piu) decays only polynomially as 
u — > 00, implies that the large deviations functional of N t /t is identically equal to 
zero between and v := lim^oo Nt/t. A key point to see this is the following: for 
a >0, 



Pr,T(y<«j = P T)T (S [at]+l > t) 

> P r , T (ti > t) ~ r 2 (3.1 



and taking the log on both sides of the inequalities, dividing by t, this implies that the 
large deviations functional vanishes at least in a = 0. Taking care of technicalities 
which involves the proof of the convexity of the functional (see [H]), one concludes. 
We go back to the computation <p k ,i' we basically identify the LHS of (J3.6P with 
the large deviations functional of the stochastic dynamics. Namely, we write: 

0/M t / 'jk(l&s 



N 2 \ N 

N 2 



E k+l (lAs) - E k (lAs),E k (lAs] 
4l • (3-9) 
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Then, using (13 .7p and ( 13. 5 p one obtains finally: 

Z(h c) 

where <p(j, r, T) is given by the RHS of (13 .7p 

2.0 



s pi 

dr / dx (fr(j(x, s'), d x e(x, s'), e(x, s')) 
Jo 



(3.10) 




Figure 3. Plot of <fi as a function of j for kt = nT 2 = 1 



4. Conclusions 



As a rigorous derivation of the relation (13. 7p for deterministic systems seems at 
the moment out of reach, it is of great importance to derive observable consequences 
which could be tested either in physical experiments or in numerical simulations. 
By analogy with the equilibrium case, we expect the lack of analyticity of <ft to 
be responsible for the occurrence of phase transitions. How those will manifest 
themselves in physical or numerical experiments is at this stage an open question. 
We note that those phase transitions are of different kind than the ones found in 
[H |2j HI [7]. In those works, they are obtained from a quadratic when optimizing 
the energy profile to get a given current. In our case, they are present in the 
system from the beginning, in the sense that <fi itself is a non-analytic function. It 
would be interesting to determine what happens to the relation giving the successive 
cumulants of the current derived in [5] which was based on a quadratic form of 0. We 
expect the relation ( 13.71) to be generic in local collisional dynamics, not only in the 
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models we described above but also in deterministic models made of tracer particles 
and fixed scatterers [121 E2 E3 1201 12] . This is because our argument is based only 
on the fact that the current is carried by a particle travelling ballistically between one 
collision and the next and on a local equilibrium hypothesis implying the occurence of 
particles with arbitrarily low speed in the system. More precisely ( 13. 7p should govern 
the macroscopic fluctuations of the system whenever it appears to be well described 
by a local equilibrium distribution. This applies to the dynamics of [101 dU US] in a 
weakly interacting regime and to the confined tracers models of [121 EH EH [201 I2T| . 
However, it does not apply to those among the models of [121 EH Ed E01 EI] in which 
the tracers are allowed to travel from cell to cell. In those models, the average over 
long times of the local transfer of energy can not be a local quantity. Indeed the local 
transfer of energy is due to particles travelling all along the system and their average 
frequency of visits to a given cell will depend on the temperatures of every scatterers 
they meet during their journey. We do not expect the behaviour described by (13. 7p 
to occur in other types of lattice Hamiltonian dynamics in which the components 
interact through a smooth potential. We emphasize also that it is not the fact that 
we replaced the deterministic collisions by stochastic ones that is responsible for the 
special behaviour of <f>. This is caused by what is left of the deterministic dynamics, 
namely the ballistic motion of the particle. In particular, one should not expect the 
same phenomenon to occur in the case of self-consistent harmonic chain [HI [91 [22] . 
The main striking feature of 0, namely its flat part, is due to the deterministic part 
that we have not modified. Therefore, one should expect this feature to be robust 
and independent of the approximations that we made for the collisional part of the 
dynamics. We note finally that a similar large deviations functional has been found 
in the context of random walks in random environments [TlTfTB] . 
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